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Интеллектуальные компьютеры — 
средство достоверного решения 
научно-технических задач 


В статье рассматриваются проблемы достоверного решения задач науки и инженерии, сводящиеся в 
конечном итоге к решению задач вычислительной математики с приближенно заданными исходными 
данными. Для решения проблем достоверности компьютерных результатов предлагается использовать 
интеллектуальное численное программное обеспечение и библиотеки произвольных разрядных вычислений, 
а также интеллектуальные параллельные компьютеры Инпарком. 


Введение 


Непрерывный рост параметров решаемых задач, постановка на компьютерах более 
полных моделей задач требуют непрерывного роста производительности компьютеров, 
который реализуется за счет распараллеливания вычислений. Эффективность решения 
задач науки и инженерии на компьютерах в большой мере зависит от достоверного ре- 
шения задач вычислительной математики, к которым в конечном итоге они сводятся. 

В настоящее время разработано огромное количество программ численного про- 
граммного обеспечения. Однако во многих случаях их применение при решении приклад- 
ных задач не приводит к получению на компьютере результатов, имеющих физический 
смысл. 

Целью данной работы является рассмотрение проблем, определяющих достовер- 
ность получаемых результатов при решении задач вычислительной математики, а также 
концепции интеллектуального компьютера, направленной на автоматизацию процессов 
исследования задачи в компьютере и ее решения с анализом достоверности получаемых 
результатов. 


Проблемы получения достоверного решения 
научно-технических задач на компьютере 


Рассмотрим проблемы достоверности компьютерных решений на одном классе задач 
вычислительной математики — решение систем линейных алгебраических уравнений 
(СЛАУ). Проблемы достоверности решений задач вычислительной математики рас- 
смотрены в [1]. 

Решение прикладных задач, как правило, начинается с создания приемлемых фи- 
зических и математических моделей. Для их построения используют различные гипотезы. 
Если эти гипотезы верны (погрешность гипотезы отсутствует или достаточно мала), то 
физическая модель правильно отражает имеющееся в прикладной задаче закономерности. 
Физическая модель может быть описана с помощью математического аппарата, напри- 
мер, некоторой системой алгебраических уравнений. 
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При описании физических моделей крайне редко возникают системы: 
АХ=Ь (0 
с точными исходными данными. 
Наиболее типичным является задание системы: 


Ах=Ь (2) 


с приближенно заданными исходными данными и указанием погрешности в исходных 
данных 


[9-4 =| 45а, 


Б-|=|АВ<в,. (3) 


Для удобства решение Х системы (1) с точными, но неизвестными исходными 
данными будем называть физическим решением, а точное решение заданной системы 
(2) — математическим решением задачи. Погрешность в решении х, вызванную неточ- 
ностью задания исходных данных, называют наследственной. Значение ее зависит как 
от погрешности исходных данных, так и от свойства матрицы системы уравнений. 

Рассмотрим некоторые примеры. 

Решение системы линейных алгебраических уравнений: 


100х, + 500х, = 1700 


: (4) 
15х, + 75,01х. = 255 
есть х! = 17, х2 = 0, 
а решение системы с приближенно заданными исходными данными: 
100х + 500х, = 1700 
(5) 


15% + 75,01х,= 255,03” 


Де: х2=3. 

При решении корректно поставленных задач, наряду с получением единственного 
классического решения, возникает необходимость в оценке наследственной погрешности 
или близости математического и физического решений. 

Верхнюю границу «относительной» наследственной погрешности точного решения 
системы можно выразить через «относительные» погрешности заданных матрицы А и 
вектора 6 следующим образом [2], [3]: 


„|. МИ“ ПАЛ, ры в 
АА М 
или 
-я АА” 
я _ ММГ АЯ ы 
Ы АМ В 
Ы 
- 5 
при условии А | | <1 и естественном предположении И <1. Обе оценки обычно 


являются мажорантными. Однако всегда можно построить пример, когда указанная 
граница достигается, т.е. оценки (6), (7) являются не улучшаемыми на всем классе не- 
вырожденных матриц. 

Из формул (6), (7) очевидно, что устойчивость решения системы к изменениям 


исходных данных в значительной степени зависит от величины с0п4 А = |“ |“ ‚ кото- 
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рая называется числом обусловленности матрицы. Если значение соп4 А невелико, то 
матрица А СЛАУ называется хорошо обусловленной. Если значение соп4 А велико, то 


матрица такой системы называется плохо обусловленной. В зависимости от способов вве- 
дения норм матрицы рассматривают несколько видов чисел обусловленности матриц [1-7]. 
Из оценок (6), (7) следует, что достоверность полученного математического реше- 
ния определяется не только обусловленностью матрицы СЛАУ, но и точностью задания 
исходных данных. 
Величину 


т = сопа А 44| [А (8) 
м 
будем называть числом обусловленности СЛАУ. Формула (8) увязывает свойства матри- 
цы системы и погрешность в задании исходных данных. В реальных задачах есть смысл 
рассматривать те системы, для которых оценки числа обусловленности т заметно меньше 
единицы, например т < 0,01. 
Повысим точность задания исходных данных в системе (4), а именно рассмотрим 
систему: 


100х, + 500х, = 1700 
15х + 75,01х,= — 255,0000037 


для которой получаем лучшее приближение к решению задачи (4) 
х! = 16,9985 х› = 0,003. 

Отметим, что при решении практических задач часто возникают СЛАУ, порядок 
которых при определенных исходных данных возрастает. Но при этом зачастую увеличи- 
вается также число обусловленности матрицы, т.е. система из хорошо обусловленной 
может превратиться в плохо обусловленную. 

Таким образом, наследственная погрешность решения определяется как значением 
величины обусловленности т, так и погрешностью задания исходных данных. И если в 
прикладных задачах нельзя повысить точность задания исходных данных, то эту задачу 
нужно переформулировать относительно других параметров и провести оценку числа об- 
условленности новой задачи. Решение СЛАУ без определения оценки числа обусловлен- 
ности может приводить к компьютерному решению, не имеющему физического смысла. 

Вторым источником погрешности решения являются ошибки округления, воз- 
никающие при вводе исходных данных задачи в компьютер (погрешность перевода 
из десятичной системы счисления в машинную, использование элементарных функций, 
которые в компьютере вычисляются по приближенным формулам и т.д.). Кроме того, 
ошибки в результатах накапливаются также при реализации вычислительных алго- 
ритмов на компьютере, поскольку все вычисления выполняются с округлением резуль- 
татов арифметических действий до некоторого фиксированного количества разрядов у 
мантиссы числа. Таким образом, в результате решения заданной системы (2) на ком- 
пьютере получают не точное решение х, а некоторое приближенное х, которое на- 
зывают компьютерным решением. 

Следуя [1], [2] будем считать, что суммарное влияние ошибок округления в пря- 
мых методах можно рассматривать как соответствующее эквивалентное возмущение 
исходных данных. 


(9) 


1 
Поэтому вычисленное машинное решение х ей 


некоторой возмущенной системы, например: 


(А+ААх® = Б+аЬ (10) 


системы (2) является точным для 
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ИЛИ (А+ЕР)х? = Б (11) 
и приближенным, не совпадающим с математическим решением системы (2). Здесь АА, 
АБ, Е — соответствующие эквивалентные возмущения, зависящие от метода решения, 
порядка системы, длины мантиссы { машинного слова и от погрешностей перевода из 
десятичной в машинную систему счисления. 
Влияние ошибок компьютерной реализации рассмотрим на следующем простом 
примере [1]. 
Пусть требуется найти решение систем линейных алгебраических уравнений: 
0135х + 0,88х› + 0,191х, -+ 0,178х, = 03516 
0,188х, + 0,262х, + 0,265х, + 0,247х, = 0,4887 
0191х, + 0,265х, + 0,281х, + 0,266х, = 0,5105. 
0,178х, + 0,247х, + 0,266х, + 0,255х, 0,4818 
С учетом того, что в математической задаче правая часть имеет четыре десятич- 
ных знака, для решения такой системы одним из прямых методов может быть интуитив- 
но выбрано одинарное машинное слово. В результате решения системы (12) методом 
[1Л-разложения (метод Холецкого) получаем машинное решение: 
0,25755808 


и _| 0,58602387 
0,63579625 | 
0,47876450 


(12) 


Вычисление невязки решения /”=р-— Ах“) с двойным машинным словом дает: 
— 0,326.10" 
0_| 09. 10° 
0,68.10—° 
0,95.10° 
С помощью фортран-программы, реализующей метод Гаусса с выбором главного 
элемента по столбцу с одинарным машинным словом, получаем машинное решение: 
0,2169916 
(_| 0,6105261 
^ | 0,6459795 |' 
0,4727254 
Если использовать С-программу метода Гаусса с выбором главного элемента 
по столбцу, то получаем следующие результаты: 
0,21508718 
\_| 061167474 | 
0,63645982 
04724402 
Однако, как нетрудно убедиться, точное математическое решение системы (12) есть 
0,4 
0,5 | 
0,6 
0,5 


Хх = 
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Решение системы (12) на компьютере с двойным машинным словом методом Гаусса 
дает компьютерное решение: 


0,4006006002160911 
0,49999998695078 
0,59999999984569 


0,5000000003221234 

почти совпадающее с точным математическим решением задачи (12). 

Из этого примера видно, что компьютерное решение задачи может сильно отличать- 
ся от ее математического решения. 

Ниже рассмотрим некоторые причины такого явления. 

В [8-19] для некоторых прямых методов приведены величины, характеризующие 
оценки норм, эквивалентного возмущения А, в частности, для метода Гаусса с выбором 
главного элемента: 


— 
А: 2 "сви, 
где { - длина мантиссы машинного слова; с! константа, независимая от и; в (п) се 
величина, зависящая от способа выбора главного элемента по всей магрице 


1 
—юп 
4 
2(п) < 1,8п тах |, |. 


Есть предположение, что для всех действительных матриц: 2(п)< птаха ‚а 


ф 


й 


Чу 


при выборе главного элемента по столбцу: 2(п)< 2"-1 пах 
5] 


Анализ суммарных ошибок округления выполненных в [1], [3], [7-9] показал, что с 
точки зрения устойчивости к ошибкам округления между прямыми методами принци- 
пиальной разницы нет. Погрешность компьютерной реализации определяется, с одной 
стороны, числом обусловленности матрицы, ее порядком, а с другой — методом реше- 
ния, способом выбора главного элемента и длиной мантиссы. Чем больше значение 1, 
тем меньше эквивалентные возмущения и тем меньше ошибки компьютерной реализации. 

Итак, при решении прикладных задач, сводящихся в конечном итоге к решению 
систем линейных алгебраических уравнений, для компьютерных моделей этих задач не- 
обходимо: 

— выявить единственное классическое или обобщенное решение; 
— определить устойчивость этого решения; 
— найти область, в которой компьютерное решение имеет физический смысл. 

Провести такие исследования могут только специалисты высокой квалификации, 
да и то не всегда. 


Интеллектуальный компьютер как средство получения 
достоверного решения 


На основе вышеизложенного можно сформулировать требования к перспективным 
интеллектуальным компьютерам [11], предназначенным для решения задач инженерии и 
науки: 

— интеллектуальный компьютер должен, без участия пользователя, исследовать 
свойства машинной модели задачи; на основании выявленных свойств, автоматически: 
строить алгоритм решения, сформировать эффективную топологию из оптимального 
количества процессоров МПМО-компьютера с точки зрения затрат машинного времени 
решения задачи, синтезировать программу параллельных вычислений на сформирован- 
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ной топологии, решить задачу и оценить достоверность решения (близость компьютер- 
ного решения к математическому решению и оценка наследственной погрешности); 

— интеллектуальный компьютер должен иметь произвольную разрядную арифме- 
тику или хотя бы библиотеку программ для работы с произвольной разрядностью для 
получения компьютерного решения, близкого к математическому решению задачи. 

На сегодняшний день численное и прикладное программное обеспечение, которое 
используется на компьютерах, возлагает оценку достоверности получаемого решения на 
пользователя, хотя это очень непростая работа, требующая глубоких знаний по вычисли- 
тельной математике, а также характерных свойств компьютерной модели задачи, мате- 
матических и технических особенностей компьютера. 

Достоверность компьютерных решений практических задач в машиностроении, 
самолетостроении и многих других отраслях, как правило, проверяется в ходе сопостав- 
ления результатов численных и натурных экспериментов, что занимает огромное время и 
увеличивает стоимость разработки объектов современной техники. 

Институт кибернетики им. В.М. Глушкова НАН Украины и Государственное научно- 
производственное предприятие «Электронмаш» совместно разработали семейство интел- 
лектуальных рабочих станций Инпарком [12-14], занимающие по своей производитель- 
ности промежуточную нишу между персональными и суперкомпьютерами, которые 
позволяют пользователям работать на суперкомпьютерах с произвольной арифмети- 
ческой разрядностью с помощью библиотек СМР [15] и МРЕК [16]. Интеллектуальное 
численное программное обеспечение исследует свойства компьютерных моделей научно- 
технических задач с приближенно заданными исходными данными и гарантирует до- 
стоверность получаемых решений. Отметим, что произвольная разрядность впервые 
аппаратно была реализована в 70-х годах ХХ в. еще на компьютерах серии МИР [17]. 


Выводы 


1. При решении научно-технических задач с приближенно заданными исходными 
данными на компьютере получаем машинные модели задач, свойства которых отличаются 
от свойств их математических моделей и, как правило, неизвестны пользователям. 

2. Необходимо исследовать характерные свойства машинных моделей задач с целью 
выбора соответствующих алгоритмов и программ решения с оценками достоверности 
результатов, а также разрядности вычислений. 

3. Интеллектуальный компьютер и интеллектуальное численное программное обес- 
печение для решения научно-технических задач обеспечивают исследование компьютер- 
ной модели задачи и, на основе полученных сведений, автоматическое построение 
алгоритма и синтез программы решения задачи с оценками достоверности на необходи- 
мой разрядности. 
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ГУ. Мосйапоъ, ГЛ. Мота, А.А. МскоаусйиЕ 

ШшеШоет Сотршег$ — а Тоо! юг пе ВеНаЫе Зо!ушто оЁ Зтепийс апа Епошеегто Ргоет5 

ТБе агасе 4еа!$ ул ргоМеплз оЁ ве тепаЫе зол оЁ ргоМептл$ ш зс1епсе ап епошееппе ата у гедисеа ю 
Фе зоУше оЁргоепл$ оР Фе сопршайопа! та®етайс$ ул арргохипаеу отуеп пу@а| даа. песет патепса! 
зой\аге ап Нбгамез оЁ ргоэтаплз Юг агФигагу ргес1юп сотршайоп$ аз уе] аз пиеШоепё рагаПе| согарщег 
шрагсот аге ргорозе4 %ю Бе изе4 Юге тезоуше оЁе геНабИИу ргоМет оЁ соппршег гези5. 


Статья поступила в редакцию 10.06.2011. 


«Штучний 1нтелект» 32011 121 


